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Abstract 

We investigate spectral functions in the vicinity of the critical 
temperature of a second-order phase transition. Since critical 
phenomena in quantum field theories are governed by classical 
dynamics, universal properties can be computed using real-time lattice 
simulations. For the example of a relativistic single-component scalar 
field theory in 2 + 1 dimensions, we compute the spectral function 
described by universal scaling functions and extract the dynamic 
critical exponent z. Together with exactly known static properties 
of this theory, we obtain a verification from first principles that the 
relativistic theory is well described by the dynamic universality class 
of relaxational models with conserved density (Model C). 



1 Introduction 



Spectral functions play an important role in our understanding of dynamic 
properties of matter in extreme conditions, ranging from heavy-ion collisions 
and physics of the early universe to laboratory experiments of ultracold 
quantum gases. In the context of collision experiments of heavy nuclei, crucial 
questions about which are the relevant dynamical degrees of freedom near 
and above the QCD transition are encoded in spectral functions. Likewise, 
transport properties such as described by shear and bulk viscosities can be 
obtained from the slope of appropriate spectral functions at small energies. 
Though these transport coefficients are computed in thermal equilibrium, 
they contain important information about time-dependent properties and 
provide essential ingredients for hydrodynamic descriptions. 

Approximate resummation techniques have been successfully implemented 
to study even far-from-equilibrium behavior and thermalization of spectral 
functions in scalar [T] and fermionic quantum field theories [2] so far. Simula- 
tions from first principles on a real-time lattice might be based on stochastic 
quantization techniques |3] , but these are still in its infancies. Lattice simula- 
tions are typically performed in Euclidean space-time, and the reconstruction 
of real-time quantities based on analytic continuation using bayesian methods 
such as MEM g] is often difficult. 

A particular challenge concerns dynamic critical phenomena near 
continuous phase transitions, such as the QCD critical point in the phase 
diagram of temperature versus baryon number density. At the transition 
point anomalously large fiuctuations signal scale invariant physics, where 
also a hydrodynamic description based on a separation of scales becomes 
inapplicable. A standard procedure for the description of dynamic critical 
phenomena is to identify the relevant hydrodynamic degrees of freedom 
away from the critical temperature with the associated dynamic universality 
class [5]. This approach successfully explains a large body of experiments in 
condensed matter physics, and the same may be expected for applications to 
relativistic heavy-ion experiments [i] though a rigorous justification remains 
an interesting problem. 

In this work we discuss the use of classical-statistical simulation 
techniques on a real-time lattice to compute the spectral function in 
the vicinity of a second-order phase transition. The universal properties 
of quantum field theories can be rigorously simulated from a classical 
counterpart in the same universality class. This is a well-established fact for 
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static properties and we explain why it also holds for dynamic properties. The 
method extends an earlier suggestion to employ the fluctuation-dissipation 
relation or Kubo-Martin-Schwinger periodicity condition |7] to compute 
spectral functions from averages of products of classical fields [8]. Non- 
equilibrium dynamics of systems near second-order phase transitions using 
classical fields were also investigated in P [ID] • As an example, we compute 
the spectral function of a relativistic single-component scalar field theory 
in two spatial dimensions. We extract universal scaling functions and the 
dynamic scaling exponent z. Together with exactly known static properties 
of this theory, we obtain a verification from first principles that the relativistic 
theory is well described by the dynamic universality class of relaxational 
models with conserved density (Model C) in the classification scheme of 
Halperin and Hohenberg [5]. 

Though only universal properties will exactly agree for the quantum 
and the classical theory, we also investigate the non-universal behavior in a 
range of temperatures around the critical temperature for which classical- 
statistical fluctuations are expected to dominate over quantum corrections. 
We observe that for temperatures of T > 1.5 the spectral function is 
approximately described by a broad Breit-Wigner form, where the non-zero 
" quasi-particle" mass and width are of the same order of magnitude. At 
around ten to twenty percent deviation from a "twin-peak" structure 
appears due to additional strong contributions from the zero-frequency 
domain of critical phenomena. At we find that the quasi-particle peak 
is suppressed and the spectral function obeys a power-law behavior at small 
frequencies and momenta. 

2 Spectral function 

For an arbitrary bosonic Heisenberg operator 0(t, x) the spectral function 
can be defined from the expectation value of the commutator as 

p(t - t', X - x', T) = ^([6(t,x),6t(t',x')]). (1) 

This expectation value in thermal equilibrium at a temperature T is obtained 
from 

(0(t,x)) = iTr(e-^/^0(t,x)) , Z = Tre-^/^ (2) 
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where the trace is taken over the Hilbert space for given Hamilton operator 
H. We also define the so-called statistical two-point function from the anti- 
commutator: 

F(t-t',x-x',T) = l({6(t,x),6t(t',x')}) -(6(t,x))(6t(t',x')). (3) 

In thermal equilibrium the spectral and statistical two-point functions are 
connected by the fiuctuation-dissipation relation or the so-called "KMS" 
condition in Fourier space 

F{u,p,T) = -i(nTicu) + ^^piu,p,T). (4) 

The Bose-Einstein distribution ?t,t(w) = l/(exp(a;/T) — 1) depends only on 
the frequency u for given temperature, and we denote the spatial momentum 
as p = |p|. In d spatial dimensions the Fourier transform reads 

f(..p.r)^/d.d^.e— .F(,x.r). (5) 

As an example, we will consider a relativistic real single-component 
scalar field theory. For the one-particle spectral and statistical functions 
0(t, x) = 0^(t, x) = 0{t,x.) with Heisenberg field operator 0{t,x.). In mean 
field type approximations the spectral function ([T]) would then be 

Poiuj,P,T) = 2msgn{u)6{uj^-p^-M\T)) , (6) 

where M{T) denotes an effective, temperature dependent mass. In the 
vicinity of the critical temperature Tc of a second-order phase transition the 
mean field behavior for the mass is given by M^(T) ~ |T — Tc|. In this case 
one observes that (|6]) obeys the scaling behavior 

po{suj,sp,s'^\T -Tc\) = s'^ po{uj,p,\T -Tc\) , (7) 

where s denotes a real, positive scale parameter. The spectral function ([6]) 
describes propagating modes with linear dispersion relation. We will study in 
detail below the emergence of slow, diffusive modes near the phase transition, 
which affect the scaling behavior significantly. 
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3 Critical dynamics 



For a discussion of the physics near the critical temperature Tc of a second- 
order phase transition it is convenient to introduce the reduced temperature 

T - T 

Tr^-^. (8) 

where we will consider T,. > if not stated otherwise. In the limit T — )■ Tc, 
for w — )■ and p — )• the spectral function can exhibit scaling behavior^ 

The critical exponent t] is the anomalous dimension, u is the correlation- 
length exponent and z denotes the dynamic scaling exponent. Comparison 
with ([7j) shows that in mean field approximations these would be 77 = 0, 
u = 1/2, and z = 1 from this quantity for the relativistic field theory. 
Dynamic universality classes can be identified, using that at sufficiently large 
distances and time scales the dynamics is described by the slowest varying 
degrees of freedom. In hydrodynamic approaches the low-energy effective 
theory includes, apart from the order parameter, the densities of conserved 
charges 0. If couplings to the latter are properly taken into account this 
leads to a mean field value 2; = 2 as will be discussed in the conclusions. 

From (Q we recover the well-known limiting cases. For instance, because 
the physics is scale- invariant, we can choose the scale parameter such that 
s^u = 1. In this case one finds from (Q for p = at T = T^. 

p(w,0,0) ~ uj-^ . (10) 
We will frequently consider ([9]) Fourier transformed with respect to time: 

p{s-%Sp,S^Tr) = S-^^-'-^'>^ pit,p,Tr). (11) 

For instance, for = this leads with s~^t = 1 to 

p(t,0,T,) = t"-^-'g(^] , (12) 



Here we consider oj > 0, t > and, due to spatial isotropy, p= |p|. 
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with the scahng function g{t/^t) = p(l) 0, {t/^tY^^'^^^)- Here we have used the 
temporal correlation length ^f, which scales with the reduced temperature as 

6 ~ r-- . (13) 

The scaling function is expected to behave as 

9 ~ . (14) 

Since the temporal correlation length f|T3|) diverges at T = Tc showing critical 
slowing down, the scaling function ( fT4l) approaches a constant and according 
to ( |T2l) the spectral function obeys a power-law behavior in time. 

For a non-zero temperature phase transition, where Tc > 0, the limit 
w — )■ implies u <^ Tc. As a consequence, in this case the fluctuation 
dissipation relation @ can be replaced by 

T T 

F{u,p,Tr) ~i — p{uj,p,Tr) . (15) 

Here the Bose-Einstein distribution leads to the classical distribution 
nxiuj) ~ T/cj for a; ^ T. Equation f|T5l) becomes exact in the limit w — )■ 
relevant for critical behavior, which leads with to the corresponding 
scaling behavior for the statistical two-point function: 

F{s'u, sp, Tr) = T(w, p, Tr) . (16) 

Equation (fT5l) can be written after a temporal Fourier transform 

p(t,p,T,) ^ -laiT(t,p,T,,). (17) 

Here the low-frequency limit for the validity of f lTSj) in the quantum theory 
translates into the long-time behavior of the correlators in the vicinity of T^.. 



4 Real-time lattice simulations 

The fluctuation-dissipation relation, given by dl]) for a bosonic quantum 
theory, relates the spectral and statistical functions in thermal equilibrium. 
It is an identity, which can be used to determine whether a theory is 
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governed by quantum or classical statistics in terms of the occupation number 
distribution nT{u}). While for a bosonic quantum theory nrp{u)) corresponds 
to the Bose-Einstein distribution, a classical-statistical bosonic field theory 
is characterized by f|T5|) with nT{uj) = T/u for all frequencies. According to 
([T5|) this classical condition is fulfilled in the quantum theory only for small 
frequencies, uj <^T . Here it is important to note that only momenta in the 
limit w — )■ 0, p — contribute to universal quantities. As a consequence, 
critical phenomena at a non-zero temperature phase transition are rigorously 
characterized by classical dynamics. 

We employ this in the following to obtain the dynamic scaling exponent 
z as well as universal scaling functions for the spectral function from 
classical-statistical lattice simulations. In a classical theory there are no 
commutation relations and the statistical two-point function corresponding 
to ([3]) becomes [11] 

- t', X - x', T) = (0(t, x)Ot(t', x'))ci - (0(t, x)),i {0{t\ x')),i . (18) 

We consider a classical Hamiltonian H{lp,7t) of the scalar field ip{t,x.) and 
its conjugate momentum field 7r(t, x). The statistical average in f|T8|) is then 
given by 

m, x)Ot(t', x'))ci = [dTTo] [dy^o] e-^/^0(t, x)Ot(t', x') , (19) 

with the classical partition function Zd = J [dvro] [dipo] e^^^^. The functional 
integral is over classical phase-space with J [dvro] [d(po] = / nxd7ro(x)d(/?o(x) 
for given fields vro(x) = 7r(t = 0, x) and v^o(x) = ip{t = 0,x.) at some arbitrary 
initial time t = 0. The time dependence of 0(t,x) = 0[ip{t,x.),7i{t,x.)] 
is then obtained from Hamilton's equations of motion for the fields ip{t,x.) 
and 7r(t, x). For known Fci(t,x, T) the classical-statistical spectral function 
Pci(i,x, T) can then be efficiently obtained from [8J 

Pei(t,P,T) = -^dtF,iit,p,T). (20) 

Alternatively, a direct computation of the spectral function in the classical 
theory would involve the more complicated computation of Poisson 
brackets [llj. Comparison of fl20|) with f|T7j) shows that near the critical 
temperature of a second-order phase transition the classical theory describes 
the long-time behavior of the corresponding quantum theory. This universal 
asymptotic behavior is characterized by critical exponents and scaling 
functions as described above. 
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5 Critical dynamics in 2+1 dimensions 



As an example, we compute the one-particle spectral function of a relativistic 
scalar field theory with one component in two spatial dimensions. It is 
described by the Hamiltonian 

H = jd\ Q + liv^r + \ -v^ (21) 

with real mass-parameter m? < and coupling A > such that A/|m| = 1 
without loss of generality for universal properties. This theory has served as 
a paradigm for the understanding of second-order phase transitions at non- 
vanishing temperature [12]. Together with exactly known static properties 
of this theory, we will use the example to obtain a verification from first 
principles that the relativistic theory is well described by the dynamic 
universality class of relaxational models with energy conservation (Model 
C) in the classification scheme of Halperin and Hohenberg |j5j . 

We extract the dynamic scaling exponent z and universal scaling functions 
from real-time lattice simulations. For this we define the field theory on a 
spatial lattice with N xN sites with periodic boundary conditions and lattice 
spacing a. Universal long-distance properties do not depend on the lattice 
spacing, but the finite volume has to be taken into account. We consider 
lattices with = 1024, 256 and 128 to verify whether convergence of the 
results is obtained for given T,,, p and u (or t). To generate thermal initial 
conditions we use a Metropolis algorithm, while the subsequent real-time 
evolution of pci{t,p, T) is obtained from Hamiltonian dynamics for the fields. 
We typically use 30-160 independently thermalized initial configurations for 
each temperature. For the Hamiltonian evolution part, we average over 
results with different initial but same relative time along the trajectory 
using time translation invariance. This increases the effective number of 
independent configurations by a factor of ~ ^t/^t-, where At is a sampling 
time of typically 1-3 lO^a. 

With 7r(i,x) = 0(t, x) the scalar field equation of motion read^ 

(pit, x) = (A - m^) ipit, x) - ^ ip\t, x). (22) 

^Replacing in (|22]) the second-order time derivative with a first-order term leads to a 
time-dependent Ginzburg-Landau equation (Model A) [5]. 
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For the numerical implementation of the equation the leap frog discretization 
is used on a cubic space-time lattice, with spatial and temporal lattice 
spacings a, to give: 

(^(t + at,x) = 2(^(t,x) - (^(t - at,x) - m2^(^(t,x) - -^(^^()f:,x) 

2 2 

+ ~l XI (^*^^' X + a«) - 2(^(t, x) + ^(t, X - cn)^ , (23) 

i=l 

where (^(t, x) = -y/ay9(t,x), = m?a?', and A = Aa is the dimensionless 
field variable, mass squared, and coupling, respectively. Here i is the unit 
vector in the spatial direction indicated by i. We solve this equation using 
m^a^ = —1 and Aa = 1 with typically at/a = 0.1. 

Since the aim is to compute dynamical behavior, we use the exactly known 
static exponents for the two-dimensional Ising universality class [12] . Taking 
the spatial correlation length exponent u = 1 and the anomalous dimension 
r] = 1/4, other static exponents are related to the two independent ones 
in terms of scaling relations such as 6 = A/t] — 1 leading to 6 = 15. We 
will verify further below as a check that our simulations reproduce well 
values of static exponents. For this, we will also couple the order parameter 
to a small external source B, adding a term Hb = — J (PxBip to the 
Hamiltonian fl2T|) . This is also used to numerically obtain an approximate 
value of the critical temperature from the vanishing of the order parameter 
(f)(T,B) = {ip{t,x.))ci = when T —t- Tc is approached from below and 
5 — )■ 0. From this we find Tc/\m\ = 4.45 ±0.05 for N = 128. To obtain more 
accurate estimates, an efficient procedure is to determine the temperature at 
which the real-time behaviour of the spectral function exhibits a power-law 
at sufficiently late times for 5 = 0, as is explained in more detail below. 
On the = 1024 lattice this gives Tj\m\ = AAU ± 0.005. If not stated 
otherwise we will consider B = and T > in the following. 

Using (1201) and (ITS]) the spectral function is given by 

p4t,P,T) = -fj d'xe-'^^^-y^ (7r(t,x)y,(0,y))ei. (24) 

Since thermal equilibrium is time translation invariant, a particular choice of 
the initial time t = in ( |2^ is irrelevant. Fig. [T] shows the spectral function 
(IMl) for different temperatures as a function of frequency u at momentum 
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Figure 1: The spectral function as a function of uj at different temperatures 
for p = 0. The approximate Breit-Wigner form at T = 2 Tc goes over to a 
"twin-peak" structure while an infrared divergence builds up as T — )■ T^. 



p = 0. For T = 2Tc its shape is approximately described by a broad Breit- 
Wigner curve whose quasi-particle peak position is at Wq/Tc ^ 0.12 while 
the width is 7o/^c = 0.036. This means that the peak appears already 
for this temperature at comparably small frequencies. As a consequence, 
for these frequency values the Bose-Einstein distribution is well described 
by riTiu)) ~ T/cj and the classical theory is expected to give an accurate 
description. Of course, results for the quantum theory can only be rigorously 
inferred in the immediate vicinity of T^. in the limit of vanishing u and p. 

As the temperature decreases, the peak position moves to smaller u. 
However, instead of a single peak continuously approaching zero near T^, 
a second peak appears before. At temperatures of about ten to twenty 
percent above Tc a twin-peak structure is visible from Fig. [H Even at only 
about one percent deviation from the critical temperature, a remnant of the 
quasi-particle peak still remains visible next to the strong peak anticipating 
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Figure 2: The critical spectral function data without connecting lines as in 
Fig. [H The double-logarithmic inset exhibits a power-law in the infrared. 



the divergence in the zero-frequency domain of critical phenomena. At 
the spectral function becomes non-analytic around u = 0: Because of the 
anti-symmetry property pc\{—uj,p,T) = —pc\(u!,p,T) the spectral function 
vanishes at w = 0, while it diverges to ±oo approaching this point from 
above and below. 

In Fig. [2] the nature of this divergence in the vicinity of becomes 
apparent from the double- logarithmic plot in the inset. The data points 
exhibit a power-law behavior towards low frequencies ranging over several 
orders of magnitude in u. Comparison of the fit value for the exponent with 
ffTU]) leads to (2 — r/)/^ ~ 0.87 and with r/ = 1/4 to the following estimate for 
the dynamic scaling exponent: 

z = 2.0 ±0.1. (25) 

Besides the statistical error, we also take into account the error in the 
determination of the critical temperature. 
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Figure 3: The spectral function as a function of time t at different 
temperatures for p = 0. Near Tc the short-time oscillatory behavior is 
followed by a power-law, before an exponential sets in for t > C,f 



The corresponding temporal evolution of the spectral function is displayed 
in Fig.|3]for different temperatures andp = 0. According to f|T2l) and f|T^ one 
expects at Tc a power-law behavior in time for diverging temporal correlation 
length C,t- If there is a small deviation from Tc then C,t stays large but 
finite and an exponential decay will always appear for sufficiently large times 
t ^ This characterization of the dynamics with respect to power-law 
and/or exponential behavior neglects possible oscillatory parts on short-time 
scales, which arise from structure of the spectral function in Fourier space 
at larger frequencies such as remnants of quasi-particle peaks seen in Figs. [T] 
and m Universal properties can only be extracted after the oscillatory short- 
time physics is damped out. Sufficiently close to Tc the three different stages 
are visible in Fig. |3l For T c:^ Tc the damped oscillations at short times 
are followed by a power-law behavior, which lasts here about two orders 
of magnitude in tTc, until at sufficiently late times an exponential decay 
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takes over. Since an infinite correlation length cannot be achieved on a finite 
computer, the power-law may be extracted before the exponential dominates. 
The data at T ~ together with 1^^ yields 1 - (2 - r/)/^ ~ 0.10 - 0.11. 
With ?7 = 1/4 this leads to an estimate for z consistent with fl25|) to two 
significant digits. Going further away from Tc the temporal correlation length 
decreases such that the exponential behavior of the scaling function g{t/C,t) 
of f lT^ quickly dominates. Accordingly, Fig. [3] exhibits after the oscillatory 
behavior an exponential decay without a clean power-law in between. 

So far, we have extracted information about the dynamic scaling 
exponent, which is obtained by putting all parameters to their critical 
values (T = p = 0) except one that is varied (w or t). Important 
universal information is also contained in scaling functions, which can be 
used to vary more than one parameter simultaneously. Here we compute 
t^-C^-^y^ p(^t,0,Tr) as a function of time t and reduced temperature T^. 
According to f lT2|) and f lT3|) . this is expected to become only a function of the 
product tTff^ sufficiently close to the critical temperature. In order to use 
this as well for an independent estimate of the value of z, first we consider 
results from different temperatures sufficiently close to Tc and compute the 
simultaneous best fit taking fit parameter for known r] and u. The 

least squares fit is based on the integral of the local width of the distribution 
of curves with different temperatures for a given tT^^. The region of this 
integration excludes very early times showing non-universal behaviour, and 
only temperatures with less than ten percent deviation from are taken 
into account. If scaling holds, then a value of z exists for which all curves at 
different are on top of each other. This method for extracting z is not very 
precise, because the quality of the fit turns out not to be very sensitive to 
changes in z. Using this measurement alone, our estimate for the dynamical 
critical exponent would he z = 2.05 ± 0.15, with a somewhat larger error 
than indicated in ( 125|) . 

The result of this procedure is shown in Fig. |H Though the curves are 
comparably noisy, which could be reduced by increasing statistics at the 
expense of more computer time, they nicely exhibit the expected scaling 
behavior for not too small tT^^ close to T^. After the non-universal short- 
time physics is damped, the scaling function exhibits an exponential behavior 
as a function of tT^^ in accordance with f lT^ . Indeed, the spectral function 
becomes well described by p{t,0,Tr) ~ exp^—t/^t) in this range. 

We use this to simultaneously fit the parameters z and C,t for each curve 
corresponding to a given temperature. The result for the temporal correlation 
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Figure 4: Shown is t^-(2-^)Ap(t, 0, T^) versus t T^^ for different temperatures. 
Sufficiently close to Tc the curves approach a scaling function, which follows 
an exponential behavior for not too early times in accordance with ( IT^ . 



length as a function of the reduced temperature is shown in Fig. |5l 
It follows a power-law, which is governed by the product of exponents uz 
according to (JT3l) . With u = 1 we find again a result agreeing with (|25i) on 
the few percent accuracy level. While all previous figures show data obtained 
on a = 1024 lattice. Fig. [5] also includes data measured on smaller lattices. 
The results indicate that we have reached the infinite volume limit within 
the given errors. 

Going away from the critical domain by sufficiently increasing T^, the 
physics becomes no longer governed by scaling behavior. To see how the 
transition to the non-universal regime happens, we have included in Fig. H] 
data for the temperature T = 1.172 Tc that corresponds to the twin-peak 
spectral function shown in Fig. [1] in Fourier space. The non-universal 
oscillatory behavior visible in Fig. H] originates from the quasi-particle 
peak of the spectral function at that temperature, while the oscillation on 
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Figure 5: The temporal correlation length as a function of the reduced 
temperature, showing a power-law behavior associated to critical slowing 
down. 



average still follow the scaling function dependence characteristic for critical 
phenomena. For even larger deviations from Tc, the non-universal dynamics 
corresponding to a Breit-Wigner like spectral function as discussed above is 
recovered. 

So far, we have considered the dynamical information which can be 
extracted from the time or frequency dependence of two-point functions in 
thermal equilibrium. A similar analysis can also be performed for higher 
n-point functions. Time translation invariance restricts all n-point functions 
to depend on relative times only. As a consequence, the one-point function 
is constant, i.e. the order parameter does not depend on time in thermal 
equilibrium. Here we consider the order parameter (j){T,B) = x))ci 
in the presence of a small external field B for temperatures both in 
the symmetric phase as well as in the phase with spontaneous symmetry 
breaking. This can be used to verify that our simulations are consistent with 
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Figure 6: Shown is 0|Tj.| " versus \Tr\ ^ for different temperatures and 
values for the i?-field using /3 = 1/8 and 6 = 15. 



the exactly known static exponents for the two-dimensional Ising universality 
class [12]. 

The scaling form of the order parameter is given by 

(f){sTr,s'^^B) = S^(j){Tr,B). (26) 

Setting sTr = sgn(Tr) one gets 

0(r„S) = \TrfMB\Tr\-^'), (27) 

where (f)±{B\Tr\~^^) = 0(sgn(Tr), -B|Tr|~'^'') denotes scaling functions above 
and below Tc, respectively. Measuring the order parameter for different 
temperatures and external field values will lead to only two curves in the 
vicinity of Tc if scaling holds. Using the known values /3 = 1/8 and 6 = 15, we 
indeed find that all data is well described by two curves with good accuracy. 
In turn, one can also use this for known exponents and unknown value of 
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Tc to obtain an estimate from fitting to the temperature where the data 
falls on two curves only. Results on a = 128 lattice are given in Fig. El 
which shows overlapping data from stepping through the temperature range 
0.88 < T/Tc < 1.1 for a set of ten temperatures with approximately equal 
step size and sufficiently small values for the 5-field. It nicely demonstrates 
scaling behavior already on rather small lattices. 

6 Conclusions 

We conclude by putting our results into the context of the standard 
classification scheme of dynamic critical phenomena, which is based on the 
identification of the relevant hydrodynamic degrees of freedom A low- 
energy description of the slow dynamics away from requires the inclusion 
of hydrodynamic modes associated with conservation laws in addition to 
the order parameter field. In a single-component real scalar field theory 
the only locally conserved quantities are energy density and momentum 
density. In particular, there is no conserved order parameter. The linearized 
hydrodynamic equations for the relativistic fluid can be reduced to two 
equations for the transverse part of the momentum density and the local 
equilibrium pressure |T3] . The hydrodynamic equation for the transverse part 
of the momentum density describes diffusive modes with quadratic dispersion 
relation. In contrast, sound wave modes at sufficiently low momenta and long 
enough time scales have a linear dispersion relation. As a consequence, they 
may be effectively integrated out since their frequencies are much higher than 
the frequencies of the dissipative modes. This means that the variations of 
pressure caused by sound excitations are fast and average out, such that 
there are no sound waves on the comparably long diffusive time scales. 
Therefore, the relevant hydrodynamic mode is the transverse component of 
the momentum density. According to the classification scheme of Halperin 
and Hohenberg [5], this low-energy effective theory should be described 
by the dynamic universality class of relaxational models with conserved 
density (Model C). For single-component theories in this universality class 
z = 2 + a/u, where a denotes the exponent for the singular part of the 
susceptibility for the conserved density. For a = in c? = 2 up to a 
logarithmic temperature dependence, one concludes: 

z = 2. (28) 
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Comparison with (1251) provides a justification from first principles that 
the relativistic theory for a single-component real scalar field is well described 
by this dynamic universality class. Some cautionary remarks are in order at 
this point. For two spatial dimensions, a value of z around two is not very 
special. For instance, Model A, describing a system with no conservation 
laws, is known to be characterized by z ~ 2.16 for d = 2 [14J. The present 
statistics is chosen to be just able to distinguish between the two. For more 
realistic applications this can be systematically improved, however, at the 
expense of more computer time. We have mentioned above that a mean 
field approximation for the critical spectral function would lead to z = 1 for 
the relativistic theory. This corresponds to propagating modes with linear 
dispersion relation. In fact, for temperatures of T > 1.5 the spectral 
function is found to be rather well approximated by a relativistic Breit- 
Wigner form with considerable width. We observe the physics of diffusive 
modes sufficiently close to the phase transition. In the hydrodynamic 
description this is a consequence of the coupling to a conserved density 
field. When the latter is taken properly into account, this consideration 
would lead io z = 2 already on a mean field level. Our simulations 
thus demonstrate that the hydrodynamic terms are properly emerging from 
the underlying microscopic model in the scaling regime near the phase 
transition. Most strikingly, for a significant intermediate temperature range 
the spectral function exhibits a twin-peak structure. The strong interplay 
between propagating relativistic modes and diffusive processes at low energies 
provides a challenge for effective descriptions in this temperature range. A 
qualitatively similar behavior would be very interesting in the context of the 
discussion of the "melting temperature" of states above the QCD transition, 
if the transition for physical strange quark mass indeed occurs close to a 
second-order phase transition, or near the QCD critical point for non-zero 
baryon number density. 
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